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ABSTRACT 

The SCARE (Structural Ceramics Analysis and Reliability Evaluation) computer 
program calculates the fast fracture reliability of monolithic ceramic compo- 
nents. The code is a post-processor to the MSC/NASTRAN general purpose finite 
element program. SCARE automatically accepts the MSC/NASTRAN output necessary 
to compute reliability. This includes element stresses, temperatures, volumes, 
and areas. The SCARE program computes two-parameter Weibull strength distribu- 
tions from input fracture data for both volume and surface flaws. The distri- 
butions can then be used to calculate the reliability of geometrically complex 
components subjected to multiaxial stress states. Several fracture criteria 
and flaw types are available for selection by the user, including out-of-plane 
crack extension theories. The theoretical basis for the reliability calcula- 
tions was proposed by Batdorf. These models combine linear elastic fracture 
mechanics (LEFM) with Weibull statistics to provide a mechanistic failure cri- 
terion. Other fracture theories included in SCARE are the normal stress aver- 
aging technique and the principle of independent action. The objective of this 
presentation is to summarize these theories, including their limitations and 
advantages, and to provide a general description of the SCARE program, along 
with example problems. 
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OUTLINE 


Designing with ceramics requires a new approach involving statistics. Inherent 
to this method is the realization that any component will have a finite fail- 
ure probability; that is, no design is fail-safe. Methods of quantifying this 
failure probability have been investigated and refined. These theories have 
been programmed into SCARE. The accuracy of the FORTRAN coding and the mathe- 
matical modeling has been verified by analytical and experimental methods. 

Using SCARE a design engineer can easily calculate the change in reliability 
due to a design change. This can lead to more efficient material utilization 
and system efficiency. 


• BASIC CONCEPTS OF BRITTLE MATERIAL DESIGN 

• STATISTICAL FRACTURE THEORIES: 

1. BASED ON STATISTICS AND OBSERVATIONS: 

UNIAXIAL WEIBULL, NORMAL STRESS AVERAGING, PIA 

2. BASED ON STATISTICS AND FRACTURE MECHANICS: BATDORF 

• THE SCARE PROGRAM 

• VERIFICATION AND APPLICATION OF SCARE 

• SUMMARY 
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CERAMICS FOR ENGINES 


Structural ceramics have been utilized for various test engine components since 
the early 1970's. This work has been sustained by the unique properties that 
ceramics offer in the areas of high-temperature strength, environmental resist- 
ance, and low density. These characteristics can result in large benefits in 
system efficiency and performance. However, the brittle nature of ceramics 
causes a high sensitivity to microscopic flaws and catastrophic fracture. The 
subsequent low reliability of ceramic components has limited their application 
in large scale engine production. 



CERAMIC PROPERTIES 

- HIGH-TEMPERATURE STRENGTH 

- ENVIRONMENTAL RESISTANCE 

- LOW DENSITY 


f 

IMPROVED EFFICIENCY AND 
PERFORMANCE 



BUT 

• HIGH CERAMIC SENSITIVITY TO FLAWS 

• BRITTLE CATASTROPHIC FAILURE 

+ 


• LOW RELIABILITY 

• LIMITED APPLICATION 
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BRITTLE MATERIAL DESIGN 


The design of ceramics differs from that of ductile metals in that ceramic 
materials are unable to redistribute high local stresses induced by inherent 
flaws. Random flaw size and orientation require a probabilistic analysis. 

The first step of a probabilistic design methodology is the determination of a 
temperature-dependent fracture strength distribution from flexural or tensile 
test specimens. From this data, the reliability for a given component geome- 
try and loading is then computed. An important characteristic of probabilis- 
tic design is that the stress distribution over the entire volume is needed. 
The design is not necessarily governed by the most highly stressed location, 
but by the entire stress field. 


• CERAMICS ARE BRITTLE AND HAVE MANY FLAWS 

• RANDOM FLAW SIZE AND ORIENTATION REQUIRE PROBABILISTIC METHOD 



SIMPLE TESTS COMPLEX PREDICTIONS 


• REQUIRES ENTIRE STRESS FIELD, NOT MAXIMUM STRESS POINT 
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STATISTICAL FRACTURE THEORIES 


A common aspect of any weakest link theory is that the component volume and/or 
surface area of a stressed material will affect its strength, whereby larger 
components result in lower average strengths. This observation led Weibull 
(1939) to propose a phenomenological model to describe the scatter in brittle 
material fracture strengths. To predict material response under multiaxial 
stresses Weibull suggested averaging the tensile normal stress in all direc- 
tions. As this approach is arbitrary and involves tedious numerical integra- 
tion, other approaches have been subsequently introduced. The most simplistic 
is the principle of independent action (PIA) model (Barnett, 1967, and 
Freudenthal , 1968). The PIA theory assumes that each tensile principal stress 
contributes to the failure probability as if no other stress were present. 
Finally, Batdorf and Crose (1974) incorporated LEFM with the weakest link 
theory to predict failure on a mechanistic basis. The model was extended to 
account for mixed-mode fracture by Batdorf and Heinisch (1978). All of these 
models are available in the SCARE program. 


I 


WEAKEST LINK SIZE STRESS STATE COMPUTATIONAL THEORETICAL 

FRACTURE MODEL EFFECT EFFECTS SIMPLICITY BASIS 


WEIBULL (1939) 

YES 

UNIAXIAL 

SIMPLE 

PHENOMENOLOGICAL 

NORMAL STRESS 
AVERAGING (1939) 

YES 

MULTIAXIAL 

COMPLEX 

PHENOMENOLOGICAL 

PRINCIPLE OF 
INDEPENDENT 
ACTION (1967) 

YES 

MULTIAXIAL 

SIMPLE 

MAXIMUM PRINCIPAL 
STRESS THEORY 

BATDORF 

(SHEAR-INSENSITIVE, 1974) 
(SHEAR-SENSITIVE, 1978) 

YES 

MULTIAXIAL 

COMPLEX 

LINEAR ELASTIC 
FRACTURE 
MECHANICS 
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WEIBULL DISTRIBUTION 


The Weibull distribution is a key ingredient for weakest link fracture theo- 
ries. As the number of flaws present in a structure is proportional to its 
volume, the failure probability increases with both the applied stress and the 
material volume. The Weibull failure probability is dependent on three statis- 
tical parameters; the Weibull modulus m, the scale parameter a Q , and the 
threshold strength a u . The Weibull modulus is indicative of strength varia- 
bility, with smaller values representing a larger variation. The scale parame- 
ter, or normalizing stress, is related to the mean strength. The threshold 
strength is usually taken as zero because of limited experimental data and the 
mathematical simplicity of the resultant equation. A similar form of the 
Weibull distribution has been developed for surface flaw induced fracture, 
with corresponding surface material parameters. 


FAILURE PROBABILITY FOR VOLUMETRIC FLAWS IN UNIAXIAL TENSION 


P, = 1 - EXP 



(a > a u ) 


P f = 0 (a < <x u ) 

WHERE 

a APPLIED TENSILE STRESS 
<x 0 SCALE PARAMETER 
<r u THRESHOLD STRENGTH 
m WEIBULL MODULUS 
V VOLUME OF STRESSED MATERIAL 
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LIMITATIONS OF THE DIRECT STATISTICAL APPROACH 


Several limitations are inherent to a purely statistical approach. One prob- 
lem occurs when the design stress is below the range of experimental data. 
Extrapolation of the Weibull distribution into this regime may yield erroneous 
results if other phenomena are present. When two flaw populations exist con- 
currently, but only one is active in the strength regime tested, the predicted 
failure probability may be low (Johnson, 1983). Further, if the threshold 
strength is not zero, the strength may be underestimated (Shih, 1980). 

Finally, an approach based only on statistics can account for stress state 
effects only in an empirical fashion. 



BASIC DATA NON-ZERO 

THRESHOLD STRENGTH 



UNDETECTED FLAW STRESS STATE EFFECTS 

POPULATIONS 
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BATDORF FRACTURE THEORY 


Recognizing that brittle fracture is governed by LEFM, Batdorf proposed that 
reliability predictions should be based on a combination of the weakest link 
theory and fracture mechanics. Conventional fracture mechanics dictates that 
both the size of the critical crack and its orientation relative to the applied 
loads determine the fracture stress. However, with ceramics the small critical 
flaw size and the large number of flaws prevent determination of the critical 
flaw, let alone its size and orientation. Instead, the combined probability 
of the critical flaw being within a certain size range and being oriented so 
that it may cause fracture is calculated. As flaw sizes correspond to strength 
levels and since strength is easier to measure than size for these microscopic 
flaws, the probability of a crack existing within a critical strength range is 
determined. This involves the derivative of the Batdorf crack density function 
which is expressed using the Weibull parameters obtained in uniaxial testing. 


P f = 1 - EXP 




dN 


d a. 


do. 


cr 


cr 


INVOLVES TWO PROBABILITIES: 

dN 

(1) P f EXISTENCE OF A CRACK IN A GIVEN CRITICAL STRENGTH RANGE] = dV da cr 

d<7 C r 

N(m,a 0 ,<r cr ) BATDORF CRACK DENSITY FUNCTION (MATERIAL PROPERTY) 
a cr REMOTE, NORMALLY APPLIED, FRACTURE STRESS OF A CRACK 


(2) P (ORIENTATION OF A CRACK SUCH THAT o e > o cr ] = — 

4ir 

a e CRACK EFFECTIVE STRESS (FUNCTION OF CRACK GEOMETRY AND ORIENTATION, STRESS 
STATE, FRACTURE CRITERION) 

0 SOLID ANGLE IN STRESS SPACE WHICH INCLUDES ALL CRACK NORMALS FOR 
WHICH a e > <r cr 
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REPRESENTATIVE SOLID ANGLES 


Fracture depends not only on the existence of a crack with a certain critical 
strength, but also on the crack orientation, the far-field stress state, and 
the crack shape. A collection of crack orientations can be described by a 
solid angle. A solid angle is measured by its subtended surface area on a 
sphere of unit radius. Therefore, the measure of the solid angle containing 
all possible crack orientations is 4*rr. The other solid angle of intest, Q, 
is defined as that angle which includes all crack normals for which a e > a cr 
(Batdorf and Crose, 1974). This assumes that fracture depends on an effective 
stress producing a singularity at the crack tip and that crack propagation 
occurs when the effective stress is greater than or equal to the critical 
stress. For uniaxial tension the effective stress is highest on cracks normal 
to the applied load. Therefore, the solid angle can be measured by the sur- 
face area of "polar caps" around the loading axis. The caps will decrease in 
size as the critical stress is increased, until a cr = o\ and Q =0. A solid 
angle representative of equibiaxial tension is also shown. 




UNIAXIAL TENSION- 
TWO CONES 


EQUAL BIAXIAL TENSION- 
EQUATORIAL BELT BETWEEN 
TWO CONES 


AS THE CRITICAL STRESS INCREASES, THE SOLID ANGLE DECREASES 
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SCARE PROCEDURAL DIAGRAM 


The previously noted fracture theories have all been implemented in SCARE 
(Gyekenyesi, 1986, and Gyekenyesi and Nemeth, 1987). The bulk of the input 
data for SCARE comes from a finite element stress analysis and, if necessary, 
heat transfer analysis of the component. MSC/NASTRA N was chosen as the pri- 
mary finite element package because of its extensive capabilities and wide- 
spread usage. The SCARE program can be readily modified to accept input from 
other methods of stress analysis. Input specifying fracture strength data 
from test specimens, usually MOR (modulus of rupture) bars, is also required. 
These specimens must be separated into two categories, those which failed 
because of surface flaws and those which fractured because of volume flaws. 

In addition, a flaw geometry and fracture criterion must be selected. The 
SCARE program then calculates the statistical parameters for both the Weibull 
and Batdorf fracture models. The survival probability is calculated for both 
volume and surface flaws for each element or subelement. Since the reliabil- 
ity of each subelement is independent, the failure probability for a given 
component under the specified load is computed. 
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AVAILABLE FLAW TYPES AND FRACTURE CRITERIA 


The diagrams below depict the flaw types and fracture criteria available in 
the SCARE program. The simple statistical fracture theories, PIA and normal 
stress averaging, do not require a specific crack geometry since they are not 
based on fracture mechanics. On the other hand, Batdorf’s fracture theory can 
be used with several different mixed-mode fracture criterion and crack geome- 
tries. The combination of a particular flaw shape and fracture criterion 
results in an effective stress equation involving far-field normal and shear 
stresses. Coplanar crack extension theories include the maximum tensile 
stress theory and the total strain energy release rate theory. In reality, a 
crack is not confined to grow within its plane, and out-of-plane crack exten- 
sion criteria are more accurate. In SCARE, these criteria are approximated by 
a simple equation (Shetty, 1987). The approach involves a semi-empirical con- 
stant which is varied to model the maximum tangential stress theory, the maxi- 
mum strain energy release rate theory, or experimental results. Because of 
the flexibility of this equation, it is the preferred model for both volume 
and surface flaws. Finally, with regard to crack geometry, semi-circular 
cracks are preferred for surface imperfections, whereas penny-shaped cracks 
best reflect the geometry of volume imperfections. 
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CODE VALIDATION 


Several sample problems have been selected from the open literature to vali- 
date the SCARE program. For example, a silicon nitride disk was spun to 
fracture (Swank and Williams, 1981). The disk geometry is shown in the first 
figure, as well as the Weibull parameters from MOR bar testing. Reliability 
calculations from the SCARE code are compared to experimental results in the 
second figure. The predictions from a shear-sensitive Batdorf fracture crite- 
rion are closest to the experimental results. However, it should be noted 
that only seven disks were tested, compared to 85 MOR specimens. This leads 
to a large degree of statistical uncertainty in the disk data and may account 
for the greater difference between experimental and predicted Pf at lower 
failure probabilities. Rigorous testing is being conducted in an in-house pro- 
gram to gather more data for code validation purposes. Furthermore, the shear- 
insensitive fracture models have been favorably tested against proprietary 
codes with these capabilities. 


ROTATING ANNULAR DISK 

VOLUME FLAW ANALYSIS 



DATA: 

NC - 132 HOT PRESSED Si 3 N 4 
m 7.65 

<7 0 74.82 MPa METER (°- 3922 ) 


N(acr) 16 ’ 30 \74 .8?MPa / PER CUBIC METER 
SHEAR-INSENSITIVE CALCULATION OF N 

fj 6.35 mm (.25 in.) 

1*0 41.275 mm (1.625 in.) 

t 3.80 mm (.15 in.) 

RPM 70 000 TO 114 000 
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PROBABILITY OF FAILURE VERSUS DISK ROTATIONAL SPEED 


PROBABILITY 

OF 

FAILURE 



CD-80-32573 


I 


3-16 


DESIGN APPLICATION 


SCARE has been used for the preliminary design of a silicon nitride mixed-flow 
rotor for application in small, high temperature engines. A single blade and 
a section of the rotor hub were analyzed using the cyclic symmetry option of 
MSC/NASTRAN. The results from the heat transfer and reliability analyses are 
shown below. Again, the shear-sensitive criterion yields a higher probability 
of failure for the same applied load. However, the regions of low reliability 
are the same for both models. 

SILICON NITRIDE MIXED-FLOW ROTOR 
TEMPERATURE DISTRIBUTION 


TEMPERATURE. 
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COMPARISON OF RISK OF RUPTURE INTENSITIES 

PER UNIT VOLUME 



PRINCIPLE OF INDEPENDENT 
ACTION, P f = 0.00082 



STRAIN ENERGY RELEASE RATE 
CRITERION, PENNY-SHAPED 
CRACK, P f = 0.0062 


2 x 10-2 
1 x 10-2 
1 x 10-4 
1 x 10-6 
1 x 10-8 
1 x 10- 10 
1 x 10-12 
1 x 10-14 
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CURRENT USAGE OF SCARE 


SCARE is a unique public-domain program which has been requested by numerous 
companies in various industries. 


INITIAL VERSION OF SCARE RELEASED TO 
(FEBRUARY 1, 1988) 


GENERAL MOTORS CORPORATION 
DETROIT DIESEL ALLISON 
FORD MOTOR COMPANY 
CHRYSLER CORPORATION 
TRW VALVE DIVISION 
NORTON-TRW CORPORATION 
EATON CORPORATION 
PDA ENGINEERING 
NASA/COSMIC SOFTWARE CENTER 
SPARTA 

WLT CORPORATION 


BOEING AEROSPACE COMPANY 
GENERAL ELECTRIC COMPANY 
BABCOCK AND WILCOX 
AVCO LYCOMING DIVISION 
DOW CHEMICAL 

TRW SPACE & TECHNOLOGY GROUP 
GARRETT TURBINE ENGINE COMPANY 
ALLISON GAS TURBINE DIVISION 
OAK RIDGE NATIONAL LABORATORY 
SUNSTRAND-TURBOMACH CORPORATION 
STRUCTURAL INTEGRITY ASSOCIATES 


INDUSTRIES INVOLVED-AUTOMOBILE, AEROSPACE, NUCLEAR, & COMPUTER SOFTWARE 
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SUMMARY 


A statistical design methodology must be used with ceramics to account for not 
only the average strength but, also, the scatter in strength. However, statis- 
tics must be supplemented with LEFM to provide a mechanistic understanding of 
the pertinent phenomenon. The improvement in failure predictions when using a 
shear-sensitive fracture mechanics based failure criterion was shown in the 
rotating disk example. The fracture mechanics/statistics combination allows 
for sound predictions when considering multiaxial stress states or concurrent 
flaw populations. This has been accomplished in the public domain finite 
element post-processor SCARE. This framework will be built on as we begin 
research on ceramic fatigue due to slow crack growth. 


• PROBABILISTIC DESIGN APPROACH MUST BE USED FOR CERAMICS 

• FRACTURE MECHANICS IS NEEDED TO ACCURATELY ACCOUNT FOR: 

- MULTIAXIAL EFFECTS 

- CONCURRENT VOLUME AND SURFACE FLAW POPULATIONS 

- SLOW CRACK GROWTH 

• SHEAR-SENSITIVE FRACTURE CRITERIA GIVE BETTER RESULTS 

• FAST FRACTURE PREDICTION CAPABILITIES ARE READILY ACCESSIBLE 
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